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SOME EFFECTS OF MASS RATIO ON THE TRANSONIC FLUTTER 
CHARACTERISTICS OF UNTAPERED 45° SWEPTBACK 
WINGS OF ASPECT RATIOS 2 AND 5-5 
By H. Neale Kelly 


SUMMARY 


A study has been made of the effects of mass ratio on the transonic 
flutter characteristics of untapered 45° sweptbaci: wings of aspect 
ratios 2 and 3-5* The experimental data, which were obtained in the 
transonic nozzle of the Langley 9- by 12- inch supersonic blowdown tun- 
nel, indicated that a given wing at a given Mach number fluttered at 
essentially a fixed dynamic pressure regardless of the individual values 
of the fluid velocity and density- In contrast to the experimental 
results, the analytical results which were based on two-dimensional 
incompressible aerodynamic coefficients indicated that, for the wings 
tested, the dynamic pressure required for flutter varied with mass ratio. 


INTRODUCTION 


Flutter theory indicates that fluid velocity and density are distinct 
variables. For certain groups of wings, however, low-speed flutter 
investigations (see, for example, refs. 1 and 2) have shown that over 
a large range of mass ratios (ratio of wing mass to fluid mass) flutter 
for a given wing occurs at essentially a constant value of the product 
of the fluid density and the square of the velocity (that is, constant 
dynamic pressure). Use of the constant -dynamic-pressure flutter concept 
greatly reduces seme of the problems associated with obtaining experi- 
mental flutter data and facilitates the correlation of such data. 


The use of the dynamic pressure as a defining flutter parameter is 
particularly advantageous at transonic Mach numbers. The variation with 
Mach number of dynamic pressure for flutter can be readily obtained in 


a variable -density wind tunnel with a single mqc 
design. 



or models of a single 
counters in flight is. 
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of course, a function only cf altitude and Mach number. Hence, by use 
of proper scale factors, flutter model results expressed in the form of 
dynamic pressure as a function of Mach number can be immediately inter- 
preted in terms of the Mach number and altitude at which flutter will 
occur in flight. If dynamic pressure cannot be used as a defining param- 
eter, however, the investigator must resort to tests of many different 
models designed to simulate various combinations of altitudes and reduced 
velocities at each Mach number. Many of the recent transonic flutter 
investigations have tacitly assumed that dynamic pressure can be used 
as a defining parameter. Available data tend to justify this assump- 
tion; however, the data are limited to a relatively small range of mass 
ratios . 

Transonic flutter studies covering a large range of mass ratios 
have been made in the transonic nozzle of the Langley 9- by 12-inch 
supersonic blowdown tunnel. The data, which were obtained in tests of 
untapered 4p° sweptback wings cf aspect ratios 2 and 3-5> sre presented 
herein as an aid in evaluating the validity of employing dynamic pres- 
sure as a basic parameter in transonic flutter investigations. 


SYMBOLS 


aspect ratio, 


2Z cos A 


b 

c 

c^ 


L h,l 


hj2 


u a 


nonaimensicnal elastic-axis position measured normal to the 
leading edge from nidchord, positive rearward, semi chords; 
speed of sound 

wing semichord normal to leading edge, ft 
wing chord normal to leading edge, inT" 
wing chord parallel to free stream, ft 

frequency of oscillation, cps 
first-bending natural frequency, cps 

second -bending natural frequency, cps 

ay 

first torsional natural frequency, c = “ cps 


GJ 


torsional stiffness about wing elastic axis, lb-ft 1 " 
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la 

l 

M 

m 


t 

V 


P 

M- 

co 


“t,l 

“fa, 2 


mass moment of inertia per unit length of wing about the 
elastic axis, slug-ft 

lengoh of wing parallel to leading edge, ft 
Mach number 

mass per unit length of wing, slugs/ft 
dynamic pressure, l/2pV^, lb/sq. ft 

no ndimens ional radius of gyration of wing about elastic axis, 
^E a /rab^, fraction semichord 

wing thickness, in. 
free-stream velocity, ft/sec 

nondimens ional wing-section center-of -gravity location 
measured normal to the leading edge from elastic axis, 
positive rearward, fraction semichord 

air density, slugs/cu ft 

mass ratio, m/jtpb^ 

frequency of oscillation, radians/sec 
first -bending natural frequency, radians/sec 


second-bending natural frequency, radians/sec 

first torsional natural frequency or a uniform cantilevered 


beam cu ajC = radians/sec 


A sweepback angle, deg 

Subscripts: 

e experimental value at flutter 

calculated value at flutter based on a three-degree-of -freedom 
analysis of form similar to that of reference 3 


R 


k 
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m measured 

c calculated using uniform cantilever team theory 


MODELS 


Plan forms of the cantilever -mounted semis pan models used in the 
present investigation are shown in figure 1. All the wings were untapered 
and swept back ^5°. One series was of aspect ratio 2; the other, of 
aspect ratio 30- 

Solid magnesium, aluminum, and steel construction were used to 
obtain the desired mass variation. In order that flutter could be 
obtained within the operating range of the tunnel, the wing stiffness 
was varied by varying the thickness (see fig. l) of the wing's hexag- 
onal airfoil section. Each wing is designated by two digits and a 
letter. The first digit denotes the aspect ratio. The letter denotes 
■che type of material (M-nagnesium, A-aluminum, S-steel) and the final 
digit denotes the various wings of a given series. 

Physical characteristics of the individual wings are listed in 
table I. The wing thickness, chord, and mass were measured; the non- 
dimensional radius of gyration was calculated by using the measured 
geometric and mass characteristics of the models. On the basis of the 
usual effective root assumptions for swept wings (see ref. 3, for example) 
and the fore-and-aft symmetry cf the hexagonal wing sections, the center 
of gravity and the elastic axis were assumed to lie along the Wing mid- 
chord line. 

Measured frequencies for the various wings are also listed in 
table I. Typical experimentally determined node lines presented in fig- 
ure 2 indicate that, for these low-aspect-ratio wings, the modes cor- 
responding to the two higher frequencies were highly coupled. (Vibra- 
tion modes of the aspect-ratio-2 wings were identified by reducing the 
span cf the higher-aspect-ratio wing in small decrements and observing 
the changes in the nodal patterns.) Nevertheless, for simplicity, the 
measured frequencies are used as uncoupled frequencies in the present 
report. Measured torsional frequencies for the aspect-ratio-2 wings 
were found to be highly erratic; therefore, calculated torsional fre- 
quencies, based on the elementary beam- theory equation for CD a given 
in the section "Symbols" and torsional stiffnesses computed from the 
data given in table I, were used in the reduction of the experimental 
flutter data for these wings (A = 2) and are listed in table I. It 
is, of course, recognized that the calculated frequencies are a very 
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poor approximation to the magnitude of the actual frequencies; however, 
they serve as a satisfactory normalizing factor in the reduction of the 
test data for the aspect-ratio-2 wing. 


TEST FACILITY 


All the experimental results reported herein were obtained in the 
transonic nozzle of the Langley 9“ "by 12-inch supersonic blowdown tun- 
nel. Top and side walls of the 7-inch-high> 10-inch-wide test section 
of the nozzle are slotted longitudinally to permit operation at and 
above sonic speed. The floor of the tunnel serves as a reflection 
plane from which semispan models are cantilever-mounted. 

Control of the tunnel stagnation pressure over the operating range 
of the nozzle (approximately 16 to 31 pounds per square inch absolute) 
is accomplished by means of a throttling valve located upstream of the 
test section. Variable and continuous regulation of the test section 
Mach number is provided by a cylindrical plunger located in the closed- 
wall part of the tunnel downstream of the test section. By extending 
the plunger into the airstream, the tunnel can be choked at any desired 
Mach number. This arrangement of throttling valve and Mach number con- 
trol plunger permits independent changes of the Mach number and stagna- 
tion pressure. In addition to the throttling valve, a quick-acting 
butterfly valve is located upstream of the test section to provide rapid 
shutdown of the tunnel. 

Condensation- free flow is assured through the use of air dryers 
and heaters which are installed upstream of the test section. Tunnel 
surveys indicate that the maximum deviation of the local Mach number 
from the test-section average varies from ±0.005 at M = 0.75 to ±0.020 
at M = 1.25. 


instrumentation 


A resistance-type, electrical strain gage was installed on the 
surface of the wing near the root (see fig. l) to establish the occur- 
rence of flutter and to indicate the frequency of the flutter oscilla- 
tion. The signal from the strain gage was amplified and fed into a 
recording oscillograph to obtain a time history of the motion. 

In addition to the strain-gage output, simultaneous and continuous 
measurements of the tunnel static pressure and the stagnation pressure 
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and temperature were recorded by the multichannel oscillograph. At flut- 
ter, independent measurements of the pressure (indicated on mercury manom- 
eters) were recorded photographically. 


TEST PROCEDURE 


The test procedure was the same as that of reference 4; flutter 
points were obtained by approaching the flutter condition by one of two 
procedures: (l) setting the Mach number control plunger in a desired 

position and varying the test-section density by increasing the tunnel 
stagnation pressure in small increments until flutter was obtained or 
(2) seating the tunnel stagnation pressure at a desired value and 
increasing (.or decreasing) the Mach number in small increments until 
the model fluttered. In many instances it was found convenient to vary 
alternately both the Mach number and stagnation pressure during the 
course of obtaining a flutter point. 


RESULTS AND DISCUSSION 


Experimental Data 


Results of the experimental investigation are listed in table II. 
(in some cases, as denoted by the blanks in the table, the flutter fre- 
quency could not be determined because of the rapid destruction of the 
model. ) Data from the table have been used in the preparation of fig- 
ure 3 which shows the variation with Mach number of the flutter-speed 

divided by the square root of the mass ratio n . The 


coefficient 


bee. 


'a 


data of figure 3, which were obtained for widely different mass ratios, 
are seen to fall into a band within which no systematic effects of mass 
ratio are evident. /it will be shown in the following paragraph that 


the parameter 


(or the companion -parameter a form used 

bco a V(I - a 

in many recent flutter reports) serves as a satisfactory correlating 
factor only if flutter at a given Mach number occurs at fixed dynamic. 

) V 

The small variation of — ■ with Mach number exhibited 

hy the data appear to "be consistent with the trends indicated hy the 
data of references 5 and 6. 
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In order to illustrate the effects of mass ratio on flutter better, 
the data have been replotted in figure 4 as the variation of the flutter- 
speed coefficient with the square root of the mass ratio for small incre- 
ments in Mach n umb er. (Effects of Mach number over the small increments 
(see fig. 3 ) are negligible.) As shewn in figure 4 the flutter-speed 
coefficients vary almost linearly with the square root of the mass ratio 
and can be approximated by straight lines which, if extended, would as 
shown by the dotted lines pass through the origins. Therefore the effect 
V_ 

of both i-i and can be expressed in terms of a single parameter 

V ^ 

— — the slope of the lines. By using the relations for h, an, 

\/ii 

r_, and q found in the section entitled "Symbols," it can be shown 

CL 

that the square of this parameter is proportional to q: 


"0 a> a \AV ^ ItGJ ' 


The terms enclosed in parentheses on the right-hand side of the equation 
are functions of the wing geometry and structure only. Therefore the 
straight lines of figure 4 indicate that, for the range of variables 
tested, flutter for a given wing at any given Mach number occurs at 
essentially a fixed dynamic press tire, regardless of the individual 
values of fluid velocity and density. 


Analytical Results 

Two-dimensional incompressible aerodynamic coefficients were used 
in an analysis similar to that of reference 3 to determine the flutter- 
speed characteristics of untapered 45° sweptback wings of aspect ratios 2, 
3.5, and 00 . in the analysis, the flutter mode was represented by the 
superposition of the first and second bending and first torsion mode 
shapes of a uniform cantilever beam. The more important parameters 
used in these calculations are of the same order of magnitude as the 
experimental values listed in table I and are as follows: 


Aspect ratio 

a or Xq, 

r 2 

HHE 

(%£/»*.) 2 


P 

2 

0 

0.25 

0.021 

0.29 

0.132 

Variable 

3.5 

0 

• 25 

.013 

. 4 o 

.422 

Variable 

00 

0 

.25 

.015 

.35 

.500 

Variable 
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The results of the calculations are presented in figure 5 in the 
form of flutter- speed coefficient as a function of the square root of 
the mass ratio. The results for the wing of infinite aspect ratio indi- 
cate a linear variation of V-p/bco a with \f\I for all values of \fil 
above about 2. Extrapolation of the linear portion of the curve for 
this wing to a zero value of \f\I indicates only a small intercept on 
the axis of V^yba> a . Thus, for this case, the analytical results indi- 
cate flutter to occur at a nearly constant value of the dynamic pres- 
sure. This result is consistent with trends shown in reference 2 for 
unswept wings of various aspect ratios. The curve for the swept wing 
of aspect ratio 3 . 5 and the lower curve for the swept wing of aspect 
ratio 2 follow the same general trend as shown by the wing of infinite 
aspect ratio at the low and intermediate values of \/\x. For the higher 
mass ratios, however, the curves for the wings of finite aspect ratio 
indicate a change in slope of the curve of the variation of Vgy'bo^ 

with v/jl. Thus, the analytical method does not predict even approximately 
a fixed dynamic-pressure flutter for the finite aspect ratio wings over 
this mass-ratio range. As an example, for the aspect-ratio-2 wing the 
analytical method predicts approximately a 55-percent decrease in q 
as the square root of the mass ratio is increased from 7-0 to 12.5 (the 
mass-ratio range covered by the experimental data for O .91 < M e < O. 98 ). 

The variation is somewhat less for the aspect-ratio-3 *5 wings. 

From a study of the flutter equations it was noted that the terms 
which produced the large decrease in the dynamic pressure required for 
flutter indicated by the calculated curves at the higher mass ratios 

involved the parameter as a multiplier. Therefore, as indicated 

by the data of figure 5 and reference 2, the large decreases in dynamic 
pressure shown for the wings of aspect ratio 2 and 3*5 would not be 
expected for swept wings of high aspect ratio or unswept wings of any 
aspect ratio. 


In order to correlate data obtained under varying test conditions, 
the results of many recent flutter investigations (see, for example, 
refs, h to 6) have been presented as the ratio of the experime ntal to 
the c alculated flutter s peed . The results of the present investigation 
indicate that cau tion should be exercised in the use and interpretation 
of this form of data presentation for low-aspect-ratlo swentback wing s 
at high mass ratios . For example, the results of the present investiga- 
tion correlated on the basis of flutter-speed ratio indicate an inde- 
pendent effect of mass ratio on the flutter-speed ratio, whereas the 
experimenta l data indlcate_that the flutter speed varies linear ly with 
the square' root of the mass ratio. 
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It is interesting to note that for the wings of aspect ratio 2 
and 3.5 at high mass ratios a second noncritical- root of the flutter 
equation is encountered. (The second root for the aspect-ratio-3 *5 
wing is encountered at flutter-speed coefficients greater than 15 and 
is not shown in fig. 4. ) Although the second root occurs at a higher 
flutter-speed coefficient, at some mass ratios it occurs at a lower 

value of the reduced velocity ( ] than that of the critical root. 

\Wr 

It was found that the second root stems from the zero airspeed torsional 
mode of vibration, whereas the critical root emanates from the second 
bending mode at the lower mass ratios and the first bending mode at the 
higher mass ratios. 


CONCLUSIONS 


The following conclusions were reached as a result of a study of 
the effects of mass ratio on the transonic flutter characteristics of 
untapered 45° sweptback wings of aspect ratios 2 and 3-5*. 

1. The experimental data indicate that, for the range of variables 
tested, a given wing at a given Mach number fluttered at essentially a 
fixed dynamic pressure regardless of the individual values of fluid 
velocity and density. 

2. In contrast to the experimental results, the analytical results 
which were based on two-dimensional incompressible aerodynamic coeffi- 
cients indicated that, for the wings tested, the dynamic pressure 
required for flutter varied with mass ratio. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 28, 1958. 
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TA3LS I-- PEYSICAL CHARACTERISTICS CP TEE XKDIVXIWAI. WUGS 


[a = 0; X,. = O-l 


slugs/ft 


0.00270 

.00225 



-* 



cps 

Cp3 


3.5- s-l 

3.5- S-2 

3.5- S-3 

3.5- S-4 

5.5- S-5 

3.5- S-6 

3.5- S-7 

3.5- S-8 

5 . 5- S-9 


3 . 5- A-l 

3.5- A-2 

3.5- A-3 

3.5- A-^ 
3-5-A-5 

3-5-A-6 

3.5- A- 


Aspect ratio 2 


0.00510 

.00^57 

.00460 

.oc5o4 

.00461 


0.248 

73 

290 

510 

.248 

75 

237 

500 

.£48 

75 

237 

500 

.247 

79 

330 

583 

.24g 

80 

242 

540 

.248 

31 

320 

582 

.248 

91 

345 

642 

.24? 

72 

276 

U97 

.248 

82 

273 

586 

.248 

75 

290 

540 

.247 

91 

298 

600 

.248 

139 

525 

975 

.250 

130 

487 

950 

.248 

124 

484 

925 

.250 

125 

500 

865 

.248 

125 

485 

900 

.248 

124 

472 

673 

.248 

122 

49S 

850 

.248 

156 

510 

9 43 

.248 

139 

525 

975 

.258 

100 

418 

900 

.247 

123 

550 

910 

.2^8 

127 

605 

920 

.249 

111 

470 

732 

.248 

145 

540 

970 

.248 

116 

565 

732 

.249 

111 1 

470 

970 

.248 

145 

540 

970 

.247 

93 

435 

773 

.247 

120 

460 

808 

.248 

129 

505 

862 

.248 

112 

508 

850 

.248 

145 

540 

925 ! 

.247 

134 

605 

S90 


Aspect ratio 3.5 


0.249 

48 

258 

4go 

.250 

44 

24p 

388 

.250 

45 

246 

389 

.248 

48 

262 

4c2 

.251 

* 5 

245 

3 88 

.250 

t4 

248 

385 

.249 

49 

259 

415 

.251 

44 

240 

584 

.248 

50 

261 

412 

.248 

65 

545 

552 

-252 | 

65 

352 

545 

.253 

67 

368 

555 

.254 

66 

3 66 

550 

.249 

67 

343 

551 

.252 

67 

358 

552 

.249 

71 

380 

570 

-253 

66 

354 

546 
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Section A - A 
(enlarged) 

Typical of wing leading 
and trailing edges 




Strain 

gage 





Tunnel floor 




Thickness ratio *010 .016 

(Streanwise) 



3. 50 

3.5-s 

3.5-A 

Steel 

Alumlnmi 

.018 

.02$ 



Figure 1.- Layout of wings. (Linear dimensions are in inches.) 
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(a) Aspect ratio 2.0. (b) Aspect ratio 5 * 5 * 

Figure 2.- Typical experimental node lines. 
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Figure 5.- Experimental variation of 


with Mach number 
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He 


(b) Aspect-ratio-3. 5 wings. 

Figure 3.- Concluded. 




(a) Aspect-ratlo-2 wings. 

Figure 4 .- Variation of the experimental flutter-speed coefficients with the square root of the 

mass ratio for small increments in Mach number. 
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